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Abstract. We investigate the total entropy production of a Brownian particle in a driven bistable system. 
This system exhibits the phenomenon of stochastic resonance. We show that in the time-periodic steady 
state, the probability density function for the total entropy production satisfies Seifert's integral and 
detailed fluctuation theorems over finite time trajectories. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.40.Jc Brow- 
nian motion - 05.70.Ln Nonequilibrium and irreversible thermodynamics - 05.40.Ca Noise 



1 Introduction 

Recent advances in non-equilibrium statistical mechanics 
address, for the first time, questions relevant to modelling 
of nanoscale machinery or to thermodynamics of small sys- 
tems [I] . Of particular interest are the so-called fluctuation 



In the present work, we probe numerically the afore 
mentioned entropy production fluctuation theorems (the 
IFT and the DFT) in the case of a Brownian particle 
placed in a double well potential and subjected to an ex- 
ternal harmonic drive. In the absence of drive, the par- 
ticle hops between the two wells with Kramer's escape 
theor^s j2^!^i::F:6:j:F:WnWWWWWWWW8\ . Theyrate rk=^ Tq-^c-^^/'^^^ [37] where tq is a characteristic 
reveal rigorous relations for properties of distribution func- 
tions of physical variables such as work, heat and en- 
tropy production for systems driven away from equilib- 
rium, where Einstein's and Onsager's relations no longer 
hold. Some of these theorems have been verified exper- 
imentally [I91[201[2T|[221IM1IM1I251[261[27^ . One of 
the fundamental laws of physics, the second law of ther- 
modynamics, states that the entropy of an isolated sys- 
tem always increases. The second law is statistical in na- 
ture. At small time scales or in small systems, one can 
observe excursions away from the typical behaviour. The 
entropy production in a given time, being a fluctuating 
quantity, can take negative values and yet the average 
entropy production over all times is always positive. En- 
tropy or entropy production is generally considered as an 
ensemble property. However, Seifert has generalized the 
concept of entropy to a single stochastic trajectory [l8l 
[3T1I32]. The total entropy production along a single tra- 
jectory involves both particle entropy and entropy pro- 
duction in the environment. The total entropy production 
is shown to obey the integral fluctuation theorem (IFT) 
for any initial condition and drive, over any finite time in- 
terval. In stationary state, a stronger fluctuation theorem, 
namely the detailed fluctuation theorem (DFT) holds. Ap- 
plications of entropy production fluctuation theorems in 
different physical contexts can be found in references [33[ 



time, AV is the energy barrier height between the two 
symmetric wells and T is the temperature of the bath. 
The random hops of the Brownian particle between the 
two wells get synchronized with the external drive if 
matches twice the frequency of the external drive. This op- 
timization condition can be achieved by tuning the noise 
intensity, and is called stochastic resonance (SR) [38 P 39 1 
HU] . Noise plays a constructive role in this case and SR 
finds applications in almost all areas of natural sciences. 
To characterize this resonance behaviour, different quanti- 
fiers have been introduced in the literature [39ll40ll4Tll42l 
I43 l l44 p 45 p 46j. The work injected into the system (or the 
thermodynamic work done on the system) per cycle char- 
acterizes SR as a bona fide resonance (40y45ii46] . Recently 
work and heat fluctuation theorems have been analyzed in 
a symmetric double well system exhibiting SR in presence 
of external subthreshold harmonic [47,48J and biharmonic 
drives. Theoretical [TTlliH) and experimental [2511^ 
studies reveal the validity of the steady state fluctuation 
theorem (SSFT) for heat and work integrated over flnite 
time intervals. In the following, we extend the study to 
fluctuation theorems for total entropy production and as- 
sociated probability density functions. 
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2 The Model 

The overdamped dynamics for the position (x) of the par- 
ticle is given by a Langevin equation [50j in a dimension- 
less form, namely 



dx 
Iti 



dU{x,t) 
dx 



(1) 



where ^(i) is the Gaussian white noise with {^{t)) — 
and ('C(i)C(i')) = 2D5{t — t'), where the noise strength 
D = ksT , Ub being the Boltzmann constant. The poten- 
tial U(x,t) can be split into two parts: a static potential 
V{x) = —\x'^ + \x^, and the potential due to external 
harmonic perturbation Vi(a;,t) = —xAsmoji. A and w 
are amplitude and frequency of the external drive, respec- 
tively. The static double well potential V{x) has a barrier 
height AV — 0.25 between two symmetrically placed wells 
(or minima) located at Xm = ±1. We have restricted our 
analysis to subthreshold forcings, Ajxml < AV . The to- 
tal potential U{x,t) = V{x) + Vi{x,t). Using the method 
of stochastic energetics (51] for a given particle trajectory 
x{t) over a finite time duration r, the physical quanti- 
ties such as injected work or thermodynamic work (W), 
change in internal energy (AU) and heat (Q) dissipated 
to the bath can be calculated. They are given by 



W 



to + T 



to 



dU{x,t) 
di 



dt. 



(2a) 



AU = U{x{ta+T),to+T)-U{x(ta),to), and (2b) 



W-AU. 



(2c) 



Equation (Pc)) is a statement of the first law of ther- 
modynamics. The particle trajectory extends from initial 
time to to final time io + '^U and Q are all stochas- 

tic quantities and we have evaluated them numerically by 
solving Langevin equation using Heun's method [47ll48l 



A change in the medium entropy {Asr 
interval t is given by 



over a time 



As 



Q 
T' 



(3) 



The nonequilibrium entropy S of the system is defined 



as 



S{t) 



dx P{x,t)lnP{x,t) = (s(r)). 



(4) 



Using Seifert defines a trajectory dependent en- 
tropy of the particle as [T51I5T] 

s{t) = -\nP{x,t), (5) 

where P{x, t) is obtained by solving the dynamical 
equation for probability density evaluated along the stochas- 
tic trajectory x{t). The change in the system entropy for 
any trajectory of duration r is given by 



As = - In 



Plixito+T),to+T) 



Poix{to),to) 



(6) 



where Po{x{to),tQ) and Pi{x{to + r),to + t) are the 
probability densities of the particle positions at initial time 
to and final time to + t respectively. Thus for a given tra- 
jectory x{t), the change in entropy As depends on the 
initial probability density and hence contains the informa- 
tion about the whole ensemble. The total entropy change 
over time duration r is given by 



Astot = Asm + As. 



(7) 



Using the above definition of total entropy production, 
Seifert has derived the IFT [HUH], i.e., 



1, 



(8) 



where angular brackets denote average over the sta- 
tistical ensemble of realizations, i.e., over the ensemble of 
finite time trajectories. 

This identity is very general and holds at any time 
and for arbitrary initial conditions. Equation ([8]) along 
with Jensen's inequality implies (Astot) > which is a 
refined interpretation of the second law. This still leaves 
open the possibility that there exists individual realiza- 
tions for which this Astot is negative. In the presence of 
external periodic perturbations, the system relaxes to a 
time-periodic steady state. In this state, a stronger de- 
tailed fluctuation theorem holds [T2 j fT8 l [3T] : 



PjAstot 

Pi-Asta 



(9) 



where Astot is evaluated over time intervals r — nTi, 
Ti being the period of the external drive. P{Astot) (or 
P{—Astot)) is the probability that the trajectory produces 
(or consumes) entropy with the magnitude As tot- 

To calculate the total entropy production, we evolve 
the Langevin system under the time-periodic force over 
many realizations of noise. Ignoring transients, we first 
find out probability density function P(x,t) in the time 
asymptotic regime. In this case, P{x, t) is a periodic func- 
tion in t with the period equal to that of the external 
drive. Having evaluated the time-periodic probability den- 
sity function, we again evolve the system trajectory. The 
heat dissipated is calculated over a period (or over a num- 
ber of periods) using ((2c|) . Thereby we obtain the change 
in the medium entropy {Asm = Q /T) . Knowing the end- 
points of each trajectory, and the time-periodic P{x,t), 
the change in system entropy As is calculated (equation 
([6])). Thus we obtain for each trajectory the total entropy 
production {Astot = Asm + As) over a single trajectory. 
To calculate the averages of the physical quantities or the 
probability distribution, Astot is obtained for more than 
10^ realizations. In the following we present the results 
where all the physical parameters are taken in dimension- 
less form. 
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3 Results and Discussions 



The work Wp calculated over a period chosen at ran- 
dom varies from realization to realization and is a random 
quantity. So are Up and Qp. However, all these quantities 
satisfy equation [2c] for each period chosen at random. The 
averages of physical quantities ((• • •)) are calculated over 
10^ realizations. 

In figure [TJ we have plotted the average work done (or 
injected work) (Wp) over a single period of the external 
drive in nonequilibrium time-periodic state, as a function 
of D for ^ = 0.1. The internal energy U being a state vari- 
able, (Up) is periodic in time and hence {^^Up) = 0. From 
equation (l2c|) we find that the average heat dissipated over 
a period (Qp) equals the average work (Wp) done over a 
period. In the same figure, average total entropy produc- 
tion over a single period, {Astot,p), as a function of D 
has also been plotted. Since entropy of system is a state 
variable, {Astot,p) = (^s„j,p) = (^). The average work 
or heat exhibits a well-known SR peak (around D — 0.12) 
consistent with the condition (at low frequency of drive) of 
matching between Kramer's rate and frequency of drive, 
which has been studied in earlier results [45ll46| . How- 
ever, peak in the {Astot,p) is not at the same D at which 
SR condition is satisfied. It is expected that at resonance, 
system will absorb maximum energy from the medium and 
being in a stationary state, will release this same energy 
back to the medium. The peak for {Astot.p) not being at 
the same temperature as that for {Qp) or {Wp) is under- 
standable as {Astot,p) = {Asm.p) = {Qp/T), i.e., peak in 
{Qp) versus T will be shifted if we plot {Qp/T) versus T. 
[Similar observations are noted in the nature of directed 
current in ratchet systems. In these periodic systems, uni- 
directional currents can be obtained in a nonequilibrium 
state in the absence of obvious bias. The average current 
exhibits a resonance peak as a function of temperature. 
Even though currents in these systems are generated at 
the expense of entropy, the peak in the total entropy pro- 
duction is not at the peak for the current. [55]] In the in- 
set of figure [ll we have plotted the relative variance of 

work ( (Wrv) = — — jjjy-^ — - — 1 and that of total entropy 

(^{Asrv) = ^^^'"ffij.lf)'' -^) ■ {Wrv) exhibits a mini- 
mum around SR condition. However, {As^v) shows a min- 
imum around the same temperature at which {Astot,p) 
shows a peak. Thus, unlike {Wrv) [30ll47ll48] . {Asrv) can- 
not be used as a quantifier of SR. This is because the 
minimmn in {Asrv) is correlated to the peak in {Astotj as 
a function of D, which itself does not occur at the value 
of the D at which resonance condition is satisfied, as dis- 
cussed earlier. It may be noted that relative variance of 
both work and total entropy production over single pe- 
riod are larger than 1, implying that these quantities are 
not self-averaging (i.e., fiuctuation dominates the mean). 
However, when the observation time for the stochastic tra- 
jectory is increased to a large number (n) of periods, the 
relative variance, which scales as n"^/^, becomes a self- 




0.1 0,2 0.3 0,4 

D 



Fig. 1. Variation of {Wp){— {Qp)) and {Astot,-p) with D, for 
A = 0.1, uj — 0.1. The curve of (Wp) versus D has been scaled 
by a factor of 8 for clarity. Inset shows the corresponding rel- 
ative variances as a function of D. 




Fig. 2. Variation of {Wp) = (Qp) and {Astot,p) with for 
A = 0.1, D = 0.15. The curve of {Wp) versus u> has been 
scaled by a factor of 6 for easy comparison. Inset shows the 
corresponding relative variances as a function of uj. 



averaging quantity, i.e., mean is larger than the dispersion 

In figure [21 we have plotted (Wp) and {Astot.p) as a 
function of uj. The injected work (Wp), exhibits a peak 
as a function of w, thus characterizing SR as a bona fide 
resonance [40 l l45 l l46] . It may be noted that the peak posi- 
tion for {Astot,p), in this case, is at the same value as that 
for (Wp) or {Qp), as expected. The inset shows the rela- 
tive variance of Astot versus frequency of external drive to 
which in turn shows a minimum at the resonance condi- 
tion. 

In figure [3l we have plotted the probability distribu- 
tion P{Astot,p) versus Astot.p, for different values of noise 
strength spanning a region of temperatures around that 
of SR {D = 0.12). For low temperature side, D = 0.06, 
P{Astot.p) exhibits a double peak structure. The peak 
around zero can be attributed to the intrawell motion. 
The small peak at higher positive values of Astot,p is 
caused by the occasional interwell transition which entails 
larger heat dissipation in the medium and contributes to 
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Fig. 3. Plots of probability distribution functions of Astot.p 
for difl'erent values of noise strength D. The fixed parameters 
are: A = 0.1, w = 0.1. 



the total entropy production via entropy produced in the 
bath, {Asm,p) = {Qp)/T. At very low temperature, D = 
0.02, the interwell motion is subdominant (particle ex- 
hibits small oscillations about the minimum). P{Astot.p) 
exhibits a single peak around {Astot,p) and the distribu- 
tion is closer to Gaussian which is not shown in the graphs. 
As temperature is increased, due to the enhancement of 
interwell motion, peak at the right increases. These mul- 
tipeaked distributions are asymmetric. The distributions 
extend to the negative side. Finite values of distributions 
in the negative side is necessary to satisfy fluctuation the- 
orems. The contribution to the negative side comes from 
the trajectories which lead to transient violations of the 
second law. For higher values of temperature, D = 0.25 
(and beyond), the peak structures merge and P{Astot,p) 
becomes closer to a Gaussian distribution. Similar ob- 
servations have been made for distributions of work and 
heat in earlier literature [47ll48j . The observed values of 
(e~'^^'°' p), from our simulations, are equal to 1.045, 1.017, 
0.980, 1.024 and 1.032, for values of temperatures D = 
0.06,0.08,0.12,0.2 and 0.25, respectively. All the values 
for (e~^^*°*''') are close to unity within our numerical ac- 
curacy, which is clearly consistent with IFT (equation ([5])). 

We have plotted P{Astot,p) and P(— Zisfot,p)e^'*'°'" on 
the same graph for two values oi D [D — 0.08 and 0.25) in 
figures [4l^a) and (b) respectively, which abides by equation 
([H]), namely the DFT. We would like to mention that the 
IFT and DFT are exact theorems for a driven Langevin 
system. Our results corresponding to figures IDJa) and (b) 
act as a check on the quality of our simulation. 

In figure [5l we have plotted probability distributions of 
changes in total entropy Astot,p, medium entropy Asrn,p 
and system entropy Asp over a single period for the pa- 
rameter values D = 0.08, w = 0.1 and A = 0.1. System 
entropy Sp (t) is a state function and its average value is a 
periodic function of time in the asymptotic regime. Thus 
average change in the system entropy over a period is zero. 
Moreover, P(Asp) is a symmetric function of Asp. The 
medium entropy is related to the heat dissipated along 
the trajectory {Asm,p = Qp/T). The nature of P{Asm.p) 




Fig. 4. (a) Both P{Astot,p) and P(-zistot,p)e"^*'°''P have been 
plotted on the same graph for D = 0.08, tj = 0.1 and ^4 = 0.1. 
These curves match to a good accuracy, thereby providing a 
cross-verification for the validity of DFT. (b) Similar plots for 
D = 0.25. Other parameters are the same as in (a). 




Fig. 

ASm 

A = 



. Plots showing distribution functions of Astot,p (green), 
(red) and of Asp (black), for D = 0.08, tj = 0.1 and 



0.1. 



is identical to that of heat distribution IS] . All these prob- 
abilities exhibit finite contribution to the negative side. 

As the observation time of the trajectory increases, 
there will be decrease in the number of trajectories for 
which Astot < 0. This is expected as we go to macro- 
scopic scale in time. To this end we have plotted in fig- 
ure [61(a) the P{Astot,np) obtained over different numbers 
(n) of cycles (or for observation times r — nTi, where 
Ti is the period of external drive). For a fixed value of 
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Fig. 6. (a) Distributions of total entropy for different numbers 
of periods (n = 1, 2, 3 and 4). The inset shows data points for 
n — 20 and the corresponding Gaussian fit with — 28.61 and 
{Astot,20p) = 14.18. Parameter values are: A = 0.1, lo — 0.1 
and D = 0.12. (b) Corresponding plots of symmetry functions 
of total entropy as a function of total entropy. 



the parameters, D = 0.12, A = 0.1 and uj = 0.1, and 
over single cycle, P{Astot,p) exhibits multi-peaked struc- 
ture which slowly disappears as we increase the period of 
observation. For larger periods, P{Astot,np) tends closer 
to being a Gaussian distribution with a non-zero positive 
mean (Astot)- We also notice that as the number of pe- 
riods increases, weight of the probability distributions to 
the negative side decreases. In the inset of figure [6l^ a) , we 
have plotted probability density of Astot taken over 20 
periods. The Gaussian fit is shown. The calculated values 
of variance, — 28.61, and of the mean, (Astot) — 14.18, 
closely satisfy the condition — 2{Astot), thereby abid- 
ing by the fluctuation-dissipation relation (see equation 
(18) of 4J). If the distribution is a Gaussian and it sat- 
isfies the DFT, then the fluctuation-dissipation theorem 

= 2{Astot) must be satisfied [4lll2U54j . The presence 
of long-time tails at large values of Astot,np are not ruled 
out (non-Gaussian nature of distribution). However, nu- 
merically it is difficult to detect them. 

In figure[6ljb), we have plotted the symmetry functions 




Fig. 7. The figure gives the plots of symmetry functions of 
medium entropy as a function of the medium entropy, for differ- 
ent numbers of periods. Parameter values are: j4 = 0.1, to — 0.1 
and D = 0.12. 



of all the curves are equal to 1, which is consistent with 
DFT. The vahdity of DFT implies IFT, but not vice versa. 
The medium entropy is extensive in time while the system 
entropy is not. Only over larger number of periods, the 
contribution to Astot,np from Asnp becomes very small 
as compared to Asm.np- This means that only over large 
time periods, Asm.np obeys a DFT relation or steady state 
fluctuation theorem as noted in the earlier literature 
48] (see figure [71). In this figure we have plotted symmetry 
functions for the medium entropy for different numbers of 
periods (n). As we increase n, the slope increases towards 
1 and hence satisfies the DFT for large n. The value of 



'In 



P(As 



P(-Astat,r.p) 



versus As 



tot,np 



for different periods. 



Irrespective of the number of periods, we find that slopes 



n over which Asm,np follows DFT depends sensitively on 
the physical parameters, unlike the DFT for Astot,np- 

SR effect is detectable even in the presence of a supra- 
threshold signal 411^, i-e., when A|a;m| > AV. Our 
studies in this regime for the total entropy production 
reveal that the distributions are broad and asymmetric. 
Moreover the multipeak structures as observed in the sub- 
threshold regime are absent. 



4 Conclusions 



In conclusion, we have studied the entropy production of a 
Brownian particle in a driven double well system which ex- 
hibits stochastic resonance. Average total entropy produc- 
tion per cycle shows a peak as a function of noise strength. 
However, it cannot be directly correlated to stochastic res- 
onance condition. Moreover, as the period of observation 
increases, contribution of negative total entropy producing 
trajectories decreases. In this nonlinear system, we have 
also analyzed these theorems in the transient regime. In 
this case, we obtain a rich structure for the probability 
distribution of trajectory dependent total entropy produc- 
tion. These results along with analytically solvable models 
will be published elsewhere [55] . 
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